Competing orders, non-linear sigma models, and topological terms in quantum 
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A number of examples have demonstrated the failure of the Landau-Ginzburg-Wilson(LGW) 
paradigm in describing the competing phases and phase transitions of two dimensional quantum 
magnets. In this paper we argue that such magnets possess field theoretic descriptions in terms of 
their slow fluctuating orders provided certain topological terms are included in the action. These 
topological terms may thus be viewed as what goes wrong within the conventional LGW thinking. 
The field theoretic descriptions we develop are possible alternates to the popular gauge theories of 
such non-LGW behavior. Examples that are studied include weakly coupled quasi-one dimensional 
spin chains, deconfined critical points in fully two dimensional magnets, and two component massless 
QED3. A prominent role is played by an anisotropic 0(4) non- linear sigma model in three space- 
time dimensions with a topological theta term. Some properties of this model are discussed. We 
suggest that similar sigma model descriptions might exist for fermionic algebraic spin liquid phases. 
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Recently there has been considerable progress in de- 
scribing the competition between various kinds of order- 
ing tendencies of quantum spin systems in two space di- 
mensions. Possible Landau-forbidden direct second or- 
der quantum transitions between Neel and valence bond 
solid phases of spin-1/2 quantum antiferromagnets have 
been described in Refs. 1. The corresponding critical 
points have been described in terms of field theories in- 
volving gapless bosonic spinon fields coupled to a fluc- 
tuating gauge field. They have hence been dubbed 'de- 
confined quantum critical points'. Other recent work has 
established the stability of 'algebraic spin liquid' phases 
of two dimensional quantum magnets^, at least within 
systematic large- iV expansions. These are critical phases 
that exhibit power law correlations in the spin'^ and other 
operators^. In the best studied cases the low energy the- 
ory is described by a conformally invariant fixed point. 
All of the existing descriptions of such algebraic spin liq- 
uids is in terms of theories of Dirac spinon fields coupled 
to a fluctuating massless gauge field. Thus these may be 
regarded as 'deconfined critical phases'. Despite the util- 
ity of the spinon-gauge description neither the spinous 
nor the gauge photon are good quasiparticles of the sys- 
tem at these deconfined critical points/phases. Indeed it 
is not clear that there is any quasiparticle description of 
the low energy fixed point theory 

A large number of open questions remain on such crit- 
ical points/phases. Is a description in terms of spinous 
and gauge fields necessary? Given the slow power law de- 
cay for various classical order parameters, is a description 
directly in terms of these orders possible? In this paper 
we study a number of related questions and suggest some 
interesting answers to these questions. 

In contrast to two dimensions competing orders in one 
dimensional spin systems are extremely well-understood: 
the classic example is the one dimensional spin-1/2 an- 
tiferromagnetic chain. These have a number of proper- 



ties reminiscent of their two dimensional counterparts. 
It is well-known that the power law phase in such spin 
chains can be described by a number of different equiva- 
lent field theories. In a semiclassical description an 0(3) 
non-linear sigma model with a topological term'"^ ob- 
tains (with 9 = tt). The 0(3) vector has the phys- 
ical interpretation of being the Neel order parameter. 
More useful is a field theory in terms of an SU{2) ma- 
trix U with a Wess-Zumino-Witten(WZW) term&. The 
SU{2) matrix U again has a simple interpretation. In- 
deed U = D + iN.a with N the Neel vector,and D the 
dimerization [i.e VBS or spin-Peierls) order parameter. 
The VBS order parameter also has the same power law 
correlations as the Neel vector in the spin-1/2 AF chain. 
Thus the WZW field theory is written precisely in terms 
of the classical order parameters that have slow correla- 
tions in this algebraic spin liquid phase. However the 
action necessarily involves the non-trivial 'topological' 
WZW term. Finally a gauge description of the spin-1/2 
chain in terms of fermionic Dirac spinous coupled to fluc- 
tuating [/(I) or SU{2) gauge flelds is also possible''-^. It 
is readily seen using bosonization techniques that this is 
exactly equivalent to the SU (2) level-1 WZW theory. 

Motivated by this we may speculate that two dimen- 
sional quantum magnets (including possible deconflned 
critical points/phases) also possess field theoretic descrip- 
tions in terms of their slow fluctuating orders provided 
certain topological terms are included in the action. Here 
we will collect a number of evidences supporting this sug- 
gestion. 

We first study quasi one dimensional systems of weakly 
coupled one dimensional spin-1/2 chains. If the inter- 
chain exchanges are unfrustrated it is expected that two 
dimensional long range Neel order will develop. Frus- 
trating interchain exchanges promote VBS ordering of 
the columnar dimer pattern with a two fold degenerate 
ground state. One approach to think about the competi- 
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tion between these two distinct ground states is in terms 
of a 2 + 1 dimensional 0(3) non-linear sigma model with 
appropriate Berry phases for the hedgehog topological 
defects^. It is known that the hedgehogs are doubled in 
a continuum description, and their proliferation leads to 
VBS order in the paramagnct^". In this approach the 
Neel order parameter is simply represented as the 0(3) 
vector. On the other hand the VBS order parameter 
has a more complicated representation and is identified 
as the hedgehog topological defect. Thus this approach 
treats the Neel and VBS orders on an unequal footing. 
Here we will show how a 'superspin' description in terms 
of a four component order parameter field may be set 
up. Three of these will simply correspond to the Neel 
vector while the fourth will be the VBS order parameter 
(which is a scalar in this case with rectangular symme- 
try). A superspin 0(4) sigma model (with some weak 
0(3) anisotropy will be derived. Interestingly this sigma 
model contains a topological term (at 6 — tt). We will 
show how the 0(4) sigma model together with this topo- 
logical term reproduces the known physics of the coupled 
spin- 1/2 chains. 

Next we consider spin- 1/2 quantum magnets on an 
isotropic two dimensional square lattice. In interesting 
recent work Tanaka and Huii have shown how a 'super- 
spin' non-linear sigma model may be derived to describe 
the Neel- VBS competition in this system. Such a sigma 
model again contains a topological term. Here we revisit 
their derivation from a slightly different perspective, and 
show the equivalence between the resulting sigma model 
theory and the spinon-gauge descriptions of Refs. Q. This 
can be accomplished in full for the case with easy plane 
anisotropy for the Neel order parameter. 

Finally we consider massless QED^ with N two- 
component Dirac fermions. For large enough N this is 
the low energy theory of the algebraic spin liquid phases 
shown recently to be stable quantum phases of SU (N) 
quantum magnets. The case iV = 4 is of direct interest 
to SU{2) magnets. Here we study N — 2 and show that 
this is again equivalent to the 0(4) sigma model with a 
topological term a.t 6 = tt. The possible duality of = 2 
QED^ to the global 0(4) model without a topological 
term was conjectured recently in Ref. 12. Our derivation 
provides partial support to this conjecture but reveals the 
presence of the topological term. 

Based on all these results we suggest that similar de- 
scriptions are possible for massless QED3 at any N. If 
correct this would be a possibly useful alternate to the 
spinon-gauge descriptions that ae currently available. 



I. ALGEBRAIC SPIN LIQUID IN ONE 
DIMENSION 

The classic example of a gapless 'algebraic' spin liquid 
is the ground state of the nearest neighbor spin-1/2 chain. 
This state has slow power law correlations for both the 
Neel and VBS correlations. In this brief section we pro- 



vide a lightning overview of its field theoretic description 
which will be useful for us later. The low energy physics 
of the spin-1/2 chain may be described by the following 
field theory 

S = So + Swzw (1) 
So = J d^x^tr{d,U^d,U) (2) 

Swzw — «r (3) 

Here U G SU{2). The number F is defined as follows. 
The field U defines a map from S*^ to S^. The volume 
in S^ bounded by the surface traced out by U defines F. 
Formally 

F = ^ y" (fye^Jktr{U-^^^UU-^^JUU-^^kU) (4) 

Here the two dimensional space-time S^ is regarded as 
the boundary of a solid ball B, and y = (yi, j/2, 2/3) are 
coordinates of B. The matrix U has been extended to 
the ball B in such a way that at the boundary it has the 
correct value for the two dimensional field theory. We 
may write U in the form 

U = (1)0 + (5) 

with ipQ + (</))^ = 1, i.e in terms of a four component unit 
vector. Then the action may be written 

So = ^{^^^)y (6) 

Swzw = J d^y^i]k^ap-f54>adi(f)pd.j4>ydk(t)S (7) 

The Neel and VBS order parameters are simply deter- 
mined in terms of the unit vector (p. We have N (f) and 
Z) - 00- The WZW field theory has global 50(4) sym- 
metry due to which the Neel vector can be rotated into 
the VBS order parameter. This symmetry which emerges 
as a property of the low energy fixed point guarantees 
that the Neel and VBS order parameters have the same 
long distance correlations. Thus the WZW field theory 
provides a 'superspin' field theory of the two dominant 
competing orders of the quantum spin-1/2 chain. Note 
that the topological WZW term is required to correctly 
reproduce the known physics of the spin chain in this 
superspin decription. 

The WZW theory must be contrasted with another 
popular field theory for the spin-1/2 chain: the 0(3) 
sigma model with a topological term. This representa- 
tion treats the Neel and VBS orders on unequal footing 
even though they eventually have the same low energy be- 
havior. Indeed the VBS order parameter is represented 
rather non-trivially as the topological charge density and 
its connections to the Neel vector are not at all obvious 
in this field theory. 

Finally it is possible to examine the spin-1/2 chain in a 
slave particle description of the spins by writing the spin 
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operator at site r as 

Sr = ^ (8) 

with fr a spinful fermionic 'spinon' field at each lattice 
site. Following standard techniques the spin Hamilto- 
nian may be recast as a U{\) (or SU{2)) gauge theory 
of these spinons coupled to a fluctuating gauge field - see 
Refs. A suitable continuum approximation can then 
be made to obtain a theory of massless Dirac spinons 
coupled to a fluctuating gauge fleld (massless QED2 in 
the U{\) case). The gauge fluctuations can then be in- 
tegrated out exactly after bosonization and the result- 
ing theory is equivalent to the WZW theory above. In 
particular as emphasized in Ref. dj the slave spinon is 
quite distinct from the true spinon known to exist in the 
spectrum of the spin-1/2 chain. Thus in this one dimen- 
sional algebraic spin liquid the spinon-gauge description 
is unnnecessary and a description in terms of the slow 
competing orders is possible. 

II. WEAKLY COUPLED SPIN CHAINS 

A. Ordered phases and transitions 

Here we consider a two dimensional system obtained as 
an inflnite array of one dimensional spin-1/2 chains cou- 
pled together by short-ranged antiferromagnetic interac- 
tions. The entire system then has rectangular symmetry 
in two dimensions. In general such coupling tends to sta- 
bilize ordering of one of the slow fluctuating order param- 
eters of the decoupled chain. For unfrustrated interchain 
couplings, coUinear Neel ordering at (tt, tt) is stabilized. 
With frustration that suppresses Neel order, columnar 
VBS ordering at (tt, 0) often gets stabilized as shown in 
recent work by Starykh and Balents--. The competition 
between the Neel state and the columnar VBS state of 
this rectangular lattice may be understood in the follow- 
ing terms. Consider the limit of reasonably strong two 
dimensional coupling. Then an I? = 2-1-1 0(3) non-linear 
sigma model description of the Neel vector fluctuations 
is useful. This model must be supplemented with ap- 
propriate Berry phases. These Berry phases only affect 
singular hedgehog configurations of the Neel vector in 
space-time- . Due to the Berry phases these hedgehog 
insertions have the same symmetry as the VBS order pa- 
rameter and hence may be identified with themiS. With 
rectangular symmetry it is known that in any putative 
continuum description the hedgehogs are doubled (note 
contrast to square lattice where they are quadrupled). 
Neel order is destroyed if the hedgehogs proliferate. Their 
Berry phases then induce VBS order. Here the resulting 
state is two-fold degenerate due to the doubling of hedge- 
hogs and may be identified with the two-fold degenerate 
columnar VBS pattern. The Neel- VBS transition in this 
case is presumably first order due to the likely relevance 
of doubled monopoles at the critical fixed point of the 



hedgehog-free theorj*i^ii^. (However there will be an in- 
teresting deconfined multicritical point). 

Now we show how these results may be reproduced in 
an alternate approach that directly attempts a 'super- 
spin' description of the two competing orders (Neel and 
VBS) in terms of a four component field. As already illus- 
trated by the analysis in d = 1, this will require incorpo- 
rating an appropriate topological term. Before supplying 
the derivation let us first understand the structure of the 
allowed superspin field theory. With a four component 
'order parameter' field, we might attempt a description 
in terms of an 0(4) non- linear sigma model with some 
0(3) anisotropy. Ignoring the anisotropy for the time 
being, the action for such a theory is defined by 

^0 = j d?xdT-^(d^$f (9) 

where is a four-component unit vector. The Neel and 
VBS order parameters are related to (j) in the same way 
as in the d = \ case. The order parameter clearly lives 
in . In three spacetime dimensions (with for instance 
boundary conditions that identify space-time with S^) 
configurations of the 4> field may be classified by an in- 
teger 'winding number' corresponding to H3(S''^) — Z. 
The winding number is expressed in terms of the <f) field 
through 

Q = j d^xdTeijktaj3-^s4>adi4>l3dj4>^dk4>5 (10) 

It is then possible to consider adding a topological Berry 
phase term to the action Sq which gives different weights 
to the different topological sectors parametrized by Q. 
Thus consider 

S^So + iOQ (11) 

The physics is clearly periodic under 6^6 + 2tt. The 
system is parity and time reversal invariant for = or 

61 = TT. 

Now we will show that a derivation of such a theory 
with a four-component field for the coupled spin chains 
naturally lead to the sigma model at 9 = n. Our deriva- 
tion will closely follow Haldane's well-known derivation 
of the 0(3) sigma model field theory with the 9 term 
for one dimensional spin chains starting from the mi- 
croscopic path integral for a single quantum spin. Here 
we will start instead from the WZW theory for each spin 
chain, and then derive a two dimensional continuum field 
theory for the field (f>. 

We start from the decoupled chain limit in which each 
chain is described by the SU{2) level 1 theory. We then 
include interchain couplings for the slow modes of this 
theory. Consider therefore the following action for the 
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coupled chains: 

S = S1 + S2 (12) 

^1 = 5]^[C/(x,r,y)] (13) 
y 

S2 = dxdT^[u${x,T,y).${x,T,y + 1) (14) 
y 

-v(j>o{x,T,y)4>a{x,T,y + 1)] (15) 

Here x is the continuous coordinate along the chain di- 
rection, y is an integer labelling the different chains, and 
U is the SU (2) matrix field which enters the description 
of each chain. The action S[U] is the SU{2) level 1 WZW 
action summed over all the chains. The term ^2 in the 
action is the interchain coupling. It has been written in 
terms of the four component vector (jj. The couplings 
M, V > 0. We have assumed there is an antiferromag- 
netic coupling of the Neel order parameter {(f)) between 
two neighboring chains, and a ferromagnetic coupling be- 
tween the VBS orders (0o)- This is completely appropri- 
ate to describe the competition between (tt, tt) Neel order 
and (tt, 0) columnar VBS order in the two dimensional 
problem. 

It is useful to first consider the simple limit where u = v 
and then to perturb around this limit. When u — v the 
action has extra global 0(4) symmetry that is broken to 
0(3) X Z2 when v. At u = w, we can write ^2 as 

^2 = -? / Y.^r{U{y + l)U{y)) (16) 

This tends to prefer a 'staggered' ordering of the U per- 
pendicular to the chains where U for even chains is the 
inverse of that for the odd chains. We therefore write 

Uiy) = g{y) y^2n (17) 
= g\y) y = 2n+l (18) 

and take the continuum limit with g assumed to vary 
slowly in both spatial directions. The first term in Si 
has a simple continuum form: 

So = ^ j cPxdTtr{d,g^d,g) (19) 

The continuum form of the WZW term requires some 
work. Substituting Eqn. 1171 we get 

^(-l)^*r[g(y)] (20) 

y 

This alternating sum can now be evaluated analogous to 
Haldane's calculation in li = 1. First write it as 

^ z{T[g{y)]-T[g{y-l)]) (21) 

y=2n 

Each term is the difference in volume in S^ bounded by 
the surfaces traced out by g{%j) and g{y — 1). With g 



smoothly varying the full sum over y can then be related 
to the total volume in S^ swept out by g{x, y, r). Specif- 
ically we get 

J dxdydT-^e,jktr{g~^digg^^djgg^^dkg) (22) 

This is precisely the 0(4) model at 6* = tt as may be 
seen by writing the SU(2) matrix g in terms of the four- 
component (j) field. 

Taking u ^ v introduces some anisotropy in the model 
which lowers the symmetry to 0(3) x Z2. Let us now 
study this model in the presence of such weak anisotropy. 
Consider defect configurations of the 0(3) vector field 
(f). In three space-time dimensions these are hedgehogs. 
In the core of the hedgehog the magnitude of 4> is sup- 
pressed. In the present model this means that the (/) 
vector points along the direction in the core. Clearly 
we can distinguish two different kinds of hedgehogs de- 
pending on whether 0o is positive or negative in the core. 
These hedgehogs are very analogous to meron-vortices in 
0(3) models with weak easy plane anisotropy. Therefore 
we will refer to these as meron- hedgehogs. Each meron 
hedgehog may be regarded as one half of a point defect 
of the 0(4) model with non-zero Q. Thus for elementary 
meron hedgehogs we have Q — ±1/2 depending on the 
sign of (/)o in the core. The topological 9 term then gives a 
Berry phase e**^ associated with each meron hedgehog. 

Let us now study various phases of the anisotropic 0(4) 
model. When the 0(3) vector orders (the Neel phase) 
the meron hedgehogs or their Berry phases play very lit- 
tle role in the low energy physics. Consider now phases 
where the Neel vector is disordered. This may be usefully 
discussed in terms of a three dimensional Coulomb gas of 
meron hedgehogs which interact with each other through 
1 /r interactions and where appropriate Berry phase fac- 
tors are associated with the hedgehogs. The action for 
this Coulomb gas takes the form 

S = S, + SB+S,nt (23) 

Sc = u^(TO^, + m%) (24) 

r 

Sb = i^'^{m+r ~ m^r) (25) 

r 

Sint = ^ mrV{r - r')mr' (26) 

rr' 

Here r is the site of some three dimensional space-time 
lattice. m±r are integers corresponding to the number 
of the two kind of ± hedgehogs at site r. The first term 
is a hedgehog core action. The second term is the Berry 
phase. In the last term rrir = m^r + rri-r is the to- 
tal hedgehog number at any site r. This term describes 
the interhedgehog three dimensional Coulomb interaction 
with V{r — r') ~ 7377- As usual an equivalent sine Gor- 
don theory is readily obtained by decoupling the inter- 
action with a potential x and summing over the integers 
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m±. It is easy to see that the resulting action takes the 
form 



S= £xK {\/xf cos(2nx) (27) 



Here e'-^ creates a strength-1 hedgehog. We thus see 
that only even strength hedgehogs appear in this sine 
Gordon theory. The Berry phases have lead to a doubling 
of the hedgehogs. The doubled hedgehogs proliferate at 
long scales in this paramagnetic phase. This leads to a 
two-fold degenerate ground state (the two ground states 
being distinguished by the expectation value of the single 
hedgehog operator). Thus phases that have unbroken 
0(3) symmetry of the anisotropic 0(4) model at 6 ^ tt 
necesarily has a two-fold ground state degeneracy. 

To understand better the nature of these 0(3) sym- 
metric phases in the particular microscopic realization of 
weakly coupled chains let us study the role of various 
discrete symmetries. In particular we must distinguish 
between parity transformations Px , Py along the chain 
direction (the x-direction) and the one perpendicular to 
the chains. Microscopically the P^. symmetry involves 
X — > —X together with — > — 0o- For the hedgehogs this 
then implies m_|_ — s- — m_, m_ — > — ™+, X ~^ ~X- Under 
Py we have y — > — y, m± —>■ —m±, x ~* ~X + tt. 

In the smooth ground states of the sine-Gordon model 
the possible values of x depend on the A„. If Ai is the 
dominant coupling then Ai < prefers e'-*^ — ±i so that 
Px is broken while Py is preserved. We may identify this 
with the (tt, 0) columnar dimer state. On the other hand 
Ai > prefers e*-^ = ±1 so that Py is broken while Px 
is preserved. This may actually be identified with the 
(0,7r) columnar dimer state [i.e vertical dimers between 
the chains). It is interesting that though we derived the 
sigma model focussing on the competition between Neel 
and the (tt, 0) VBS state it is still capable of describing 
the (0, tt) VBS state as well. 

All of this is completely consistent with results known 
from the earlier direct analysis of the 0(3) sigma model 
with Berry phases appropriate for a two dimensional lat- 
tice with rectangular symmetry. Thus the superspin 0(4) 
formulation with the Q term correctly captures the Neel- 
VBS competition in this system. 

Finally we note that spin- 1/2 magnets on lattices with 
rectangular symmetry are not expected to have trivial 
featurless paramagnetic phases that also preserve all lat- 
tice symmetries. Non-trivial paramagnetic phases such 
as spin liquids with topological order are however possi- 
ble. This then implies that the 0(4) model at 6' = tt in 
the presence of 0(3) x Zi anisotropy will also not possess 
simple phases which break no symmetry. Any such sym- 
metry unbroken phase must necesarrily also have some 
hidden order (such as topological order). 



III. DECONFINED CRITICALITY IN THE 
TWO DIMENSIONAL SQUARE LATTICE 

We now turn our attention to spin-1/2 quantum mag- 
nets on isotropic two dimensional square lattices. Here 
apart from the well-known Neel state, paramagnetic VBS 
states are again possible. In contrast to the quasi-two di- 
mensional case here the VBS state will have a four fold 
ground state degeneracy. A superspin description of the 
Neel- VBS competition was recently derived by Tanaka 
and Hu^^ . Here we first very briefiy review their deriva- 
tion from a slightly different perspective. In the pres- 
ence of some easy plane anisotropy for the Neel vec- 
tor we explicitly show the equivalence to the spinon- 
gauge field theory proposed in Ref. 'J for the decon- 
fined critical point between these two ordered phases. 
The appropriate spinon-gauge field theory involves gap- 
less bosonic spinons coupled to a gapless non-compact 
U{t) gauge field and has been dubbed the non-compact 
CPi (NCCPi) modeli^. 

A. Tanaka-Hu superspin field tlieory 

Consider a half-filled Hubbard model on a two dimen- 
sional square lattice with tt fiux through each square pla- 
quette described by the Hamiltonian 

= - ^ (4cr' + /i-c) + C/^nr(nr - 1) (28) 

<rr'> r 

Here Cra are spinful electron operators at sites i of a 
square lattice. The t„' are chosen to have tt flux through 
each plaquette. A specific choice is trr' = — (— l)^t for 
horizontal bonds and trr' = t for vertical bonds. The U 
term is the usual onsite Hubbard repulsion. At [/ = the 
band structure is well-known and consists of four Fermi 
points (at (±-|,±^)) in the Brillouin zone of the square 
lattice. The electron dispersion is linear in the vicinity of 
these points. The low energy physics may be described 
in terms of a continuum Dirac theory that focuses on the 
modes near these Fermi points. For r = {x, 2y) write 

CP = e^^i->i2+e^^^-">2i (29) 
cp+y = *(e*'^^->ii+e*'^->22) (30) 

with Ki = (f , f ), and K2 = (-f , f ). Then the contin- 
uum Hamiltonian reads 

J d^xi)^ i-idxT' - idyT-") ^ (31) 

Here ip = ipaa and each ipaa is a two-component Dirac 
spinor. The Pauli matrices r act on the Dirac index. The 
index a = 1 , 2 labels the node and a labels the physical 
spin. The corresponding action (after letting y — > —y) is 

S = j (fxdr^P {-iTydr - idxT'= - idyT^) (32) 
= j Sxi){-iTidi)'<p (33) 
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where we have defined 'ij) = itp^Ty. 

At small U the interactions renormalize to zero and 
a semi-metallic state is stable. As U is increased there 
will be a metal-insulator transition where the fermions 
acquire a charge gap. In the resulting insulator it is 
expected that in some parameter ranges a simple Neel 
state will emerge. In other parameter ranges a valence 
bond solid state may well be stabilized. Either insulator 
may conveniently be described starting from the Hub- 
bard model in a Hartree-Fock description that builds in 
a non-zero value for the appropriate order parameter. To 
describe the competition between these two kinds of or- 
dered insulators, we will introduce fields that couple to 
either order parameter and integrate out the fermions to 
obtain effective sigma model descriptions. It is easy to 
see that the Neel vector coresponds to the fermion bi- 
linear iil)anyil) where are Pauli matrices in the node 
index. The VBS order parameter (vx,Vy) corresponds 
to Vx — iip^ix^'j = j'/'/^zV'- To access either insulating 
phase we therefore include the terms 



UN (ipdnyijjY + (ijjuxil^Y + {'^^^zi>y 



(34) 



To begin with we will assume that un = Uy = u for 
simplicity. Later we will add anisotropy between the 
Neel and VBS fluctuations. This term may be decou- 
pled with a 5-component Hubbard Stratonovich field 
(j) = {N, <j)vy,4>vx)- We therefore get 



S 



'iTidi + icjy.T ) ip + 



2u 



(35) 



Here F are a set of five 4x4 matrices defined through 



Fl,2,3 
F4 

F., 



fJ-z 



(36) 
(37) 
(38) 



Note that [Fa,F^]+ = 26ap and Fg = -F1F2F3F4. 
Equivalently the phase diagram and universal aspects of 
the various phases will be preserved by restricting the 
(j) vector to have unit magnitude. We thus consider the 
action 



S 



(39) 



The Hartree-Fock aproximation to the insulating phases 
consists of giving (j) some non-zero mean value which gaps 
out the fermions. To access these insulators while includ- 
ing fluctuation effects we integrate out the '0 assuming 
that m is large. The resulting fermion determinant has 
been evaluated by Abanov and Wiegmann!'^^ and gives 



S[cf>] 



(40) 



This takes the form of a non-linear sigma model for the 
superspin vector that combines the Neel and VBS or- 
der parameters. The crucial feature is the presence of 



the second term - this is a Wess-Zumino-Witten term 
for the 5-component unit vector field cj) and is defined as 
follows. The field (j) defines a map from to 5"*. The 
fraction of the total volume of that is bounded by 
the hypersurface traced out by </) defines F. Specifically 
let (j){x,u) be any smooth extension of (j){x) such that 
i{x, 0) = (1, 0, 0, 0, 0) and ^{x, 1) = ^{x). Then 

F = / du [ d^xeaf:ijSK4'adx4>i3dy(t)jdT4>sdu4>K 



87r2 

(41) 

Thus the Neel- VBS competition for spin-1/2 magnets on 
an isotropic two dimensional square lattice is described 
by this 5*0(5) superspin non-linear sigma model with the 
extra topological WZW term. This sigma model must be 
supplemented with an anisotropy that reduces the sym- 
metry to 50(3) X U{1) where the 50(3) corresponds to 
spin rotation and the U{1) rotates between the two com- 
ponents of the VBS order parameter. Further microscop- 
ically this U{1) symmetry must have further anisotropy 
that reduces the symmetry to Z4. 

In their original derivation Tanaka and Hu viewed the 
TT flux state as a ferniionic spinon mean field theory for 
Heisenberg spin magnets. Their subsequent calculations 
are identical to what we described above. However from 
that point of view it becomes important to include gauge 
fluctuations before establishing a firm connection to the 
original quantum magnet^''. Here we have sidestepped 
this issue by considering the Hubbard model on the 
square lattice with tt flux as our starting point. There is 
no question of introducing extra fluctuating gauge flelds 
from this point of view. Note that once the system en- 
ters the insulating phase the spin physics at energy scales 
below the charge gap may be described (in principle) in 
terms of a short ranged spin model of spin-1/2 moments. 
Thus the 7r-flux Hubbard model as a microscopic starting 
point is very convenient to derive an effective theory to 
describe the Neel- VBS competition in two dimensional 
quantum magnets. 



B. Equivalence to NCCP^ field theory 

We now show the equivalence of the superspin sigma 
model with the WZW term to the NCCF^ model in Ref. 
for the deconfined critical point between the Neel and 
VBS states. We will do this in the specific case where 
there is large easy plane anisotropy on the Neel vector so 
that it points primarily in the xy plane in spin space. To 
handle this case we write 



(0,H) 



(42) 



with H — iNj_, (f)yy, (j)yx) ^ four component unit vector. 
Here N±_ is a two component vector which denotes the 
direction of the Neel vector in spin space. With this 
restriction the Wess-Zumino term in the action can be 
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evaluated straightforwardly and simply becomes 

Swzw = -ii^Q (43) 

where Q is the by now familiar index describing the wind- 
ing of the unit four vector field 11 in three space-time di- 
mensions. Thus we appear to be back to the 0(4) sigma 
model ai — IT but now with 0(2) x 0(2) anisotropy for 
the n field^i 

The topological term again plays a crucial role. To ex- 
pose its role let us warm up by considering some simple 
configurations of the fl field. It will be useful to view it 
as a two component complex vector z = (zi, Z2) with zi^2 
complex numbers satisfying 1^1 P + = 1. The topo- 
logical index Q is readily writen in terms of derivatives 
of z. An equivalent form is obtained by considering the 
"vector potential" 

ai = -iz^diZ. (44) 
Some algebra shows that 

(5=^ J d^xe^jkaidjak (45) 

As a topological index Q is invariant under smooth con- 
tinuous deformations of z(x). As a particular example a 
smooth gauge transformation z —> e^^^^-'z changes at to 
Ui + diO. The integral for Q above is clearly invariant 
under such a gauge transformation. 

Consider first configurations where zi — |zi|e*^i'^^) 
with 61 a smooth single valued function of x, and Z2 
arbitrary. This means that that there are no vortices 
anywhere in the complex field zi . Then the gauge trans- 
formation z — >• e~**i'-^-'2; makes zi real everywhere. It is 
clear then that for such configurations Q — [as, with 
zi real one component of 11 vanishes everywhere). The 
same result obviously also holds for configurations where 
there are no vortices anywhere in Z2 but zi is arbitrary. 
Thus Q is non-zero only if there are vortices in both zi 
and Z2- 

To construct a configuration with non-zero Q consider 
the map 

n(r, 0, (/>) = [cos{a{r)), sin{a{r)er) (46) 

with (r, 0, p) the spherical coordinates for space-time and 
a{r) some smooth function of r satisfying a{Q) — 
and a (00) = tt. Here eV is a radial unit vector de- 
fined by {cos9, sinOcoscf), sin9sin(j)) (the first component 
is along the r direction). This configuration clearly has 
(5 = 1 and is a simple generalization of the familiar 0(3) 
skyrmion in two dimensions. Let us identify the location 
of the vortex lines of zi 2 in this configuration. When 
6* = 0, z = (e*",0). Similarly when 6* = tt, z = (e-*",0). 
Thus along the line A (see Fig. which runs along the 
T axis at the spatial origin z — (e'^^ , 0) with 9i varying 
from — TT to TT. Similarly at r such that a{r) — ^ and 
= ^ {i.e along the curve B) we have z = (0, e*"^). Thus 
we may identify A with a vortex line in Z2 and B with a 



A 




FIG. 1: Configuration with Q = 1 showing the location of the 
vortices in 21 and Z2. The curve A is the vortex in Z2 and B 
is the vortex in zi. 



vortex line in zi . In this configuration these two oriented 
vortex lines have a non-trivial linking number 1. 

These considerations illustrate a general mathematical 
result which we will henceforth use: The topological in- 
dex Q may be identified with the linking number between 
oriented vortex loops in zi and those in Z2. For the 9 = tt 
model this implies that whenever a zi-vortex moves fully 
around a Z2-vortex (or vice versa) there is a phase factor 
of TT. In other words the zi and Z2 vortices are mutually 
non-local and see each other as sources of 7r-fiux. 

It is now straightforward to see the equivalence to the 
easy-plane NCCP^ action. Indeed the vortices in zi are 
just the 27r vortices in the easy plane Neel order param- 
eter. The TT phase shift that these vortices acquire when 
taken around the vortices in Z2 suggests that the latter 
may be thought of as spin-1/2 spinous. This is indeed 
consistent with the expected structure of the vortices in 
the VBS order parameter—. To explicitly demonstrate 
the equivalence to the NCCP"'^ action let us go to a dual 
representation in terms of the vortices in zi and Z2. A 
lattice action which incorporates the non-trivial mutual 
statistics of these two vortices is readily written down: 

S ^ Sn + St2 + Sa + Scs (47) 
Sti = X! '^y'Cos(V(/)i - ai) (48) 

'5(2 = X! M/JCOs(V02 - a2) (49) 
<IJ> 

Sa = (v X ai)' + (v X 02)' (50) 

p 

Scs ^ T(l-A*/./)(l-n^) (SI) 

</./> p 

Here are the sites of a cubic space-time lattice 

and I, J,... the sites of the corresponding dual lattice. 
The fields (j)i^2 are the phases of the vortices in zi and 
Z2 respectively. These vortices are coupled minimally to 
their respective non-compact U{1) gauge fields ai_2- The 
Ciji^/j are Z2 gauge variables that take values ±1. The 
last term is an Ising mutual Chern-Simmons terrain that 
incorporates the mutual tt statistics of the two vortices. 
This will be made explicit below. The symbol Yip rep- 
resents a product of (T^ over the four links of the space- 
time plaquette pierced by the dual link < IJ >. 
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To check that the Scs term correctly implements the 
mutual statistics of the two vortices it is convenient to go 
to a Villain representation. For instance for the 2-vortex 
the hopping term becomes 



Sr. 



J2 "(-^2 



<IJ> 



d2 + ^il-fiij) (52) 



A sum over the integer vortex currents j2 is implied. The 
integration over 02 as usual implies the current conser- 
vation condition 



V.J2 = 



(53) 



The sum over the Ising variable /i/j may also be per- 
formed and leads to the constraint 



(-ir=n^ 



(54) 



Thus the presence of a 2-vortex is seen as tt flux by the 

1- vortex. Thus the term Scs precisely implements the 
mutual TT statistics of the two vortices. 

In the Neel phase the vortex e"^"^ is gapped while the 

2- vortex is condensed (with the reverse being true in the 
VBS phase). At a putative second order transition be- 
tween the two, we might expect that both vortices stay 
critical. A standard duality transformation on one of the 
two vortex degrees of freedom enables us to demonstrate 
that this action is indeed in the easy plane NCCP^ uni- 
versality class near this transition. This is elaborated in 
Appendix ^ 

Note that the easy plane deconfined critical fixed point 
also has 0(2) x 0(2) symmetry^ consistent with the su- 
perspin field theory discussed in this section. One of these 
0(2) symmetries is realized simply as a uniform phase ro- 
tation of both the matter fields e'^± of Appendix 1X1 The 
other 0(2) is realized non-trivial as a topological sym- 
metry associated with the conservation of the photon. 



IV. SIGMA MODEL DESCRIPTION OF 
MASSLESS QEDs 

In this section we will consider massless QED^ in D — 
2-1-1 dimensions. This consists of A'^ two-component 
Dirac fermions ip coupled to a fluctuating non-compact 
C/(l) gauge field a with the action 



S 



1 



d xi) {n{-idi - ai)) ■0 + TT^ {UjkdjakY (55) 



The Ti are Pauli matrices. This action has a global 
SU{N) symmetry asociated with unitary rotations of V' 
and a hidden global U{1) symmery associated with gauge 
flux conservation. For A = 4 this is the low energy de- 
scription of the dRVB algebraic spin hquid of SU (2) in- 
variant spin-1/2 models that has been much discussed re- 
cently. Neither the Dirac fermions nor the gauge photon 
are well-defined quasiparticle excitations of the system 



when its low energy properties are controlled by a scale 
invariant fixed point (as happens generically for large 
enough A^ and at special multicritical points for small 
A^). For any N > 2 this scale invariant theory has slow 
power law correlations for a set of gauge invariant fermion 
bilincars that transform as adjoints under the SU{N). 
Can we dispense with the gauge description completely 
in favor of some sort of sigma model decription in terms 
of these slow fluctuations? As discussed in previous sec- 
tions this is certainly possible in various other analogous 
problems but required inclusion of a topological term in 
the sigma model action. Here we study the specific case 
A^ = 2inD = 2-|-l and show that a sigma model de- 
scription is again possible provided topological terms are 
included. 

For N — 2 the fermions transform with spin-1/2 under 
the global SU{2) symmetry. It is expected that the gauge 
invariant SU (2) vector ipaip will have slow correlations, 
and perhaps might even generically order. To expose its 
effects let us add the following term to the action 



S4 



(56) 



with u > 0. This term may be decoupled with a fluctu- 
ating A' field to get a term of the form 



{Nf 
2u 



iN.Tpaip 



(57) 



Equivalently the phase diagram and universal aspects of 
the various phases will be preserved by restricting the 
A vector to have unit magnitude. We thus consider the 
action 



Sf + Sa 



(58) 

Sf ^ J d'^xtp {Ti{—idi — at) + imfi.a) ip (59) 

(60) 



1 2 



with m > 0. The fields V'l o-, fi are all to be integrated 
over. We begin by first doing the quadratic integral over 
the fermion fields. We assume space-time to be S^ . The 
resulting fermion determinant has been calculated (in a 
1/771 expansion) by Abanov and Wiegmanni^ and leads 
to the following remarkable result: 



S ^ So + 
So ^ I d^x- {d.h) 



Sj + Sh + Sa 

9 



5, 



d xiJiGi 



Sh = —iTTH[ 



(61) 
(62) 

(63) 
(64) 



Here g ~ 1/m in Sq. The Ji is the topological current 
density of the h field, i.e it is the skyrmion current den- 
sity. Formally 



Ji = -!-e„fe77.9o77 X dkfi 
air 



(65) 
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The last term of the action involves an interesting topo- 
logical invariant H corresponding to 113(5'^) = Z and is 
known as the Hopf term. H is an integer that distin- 
guishes different space-time configurations of the h field. 
H is conveniently written in terms of a CP^ spinor z 
associated with the field n (defined through n — z'^az). 
From z form the SU{2) matrix U — [z ~ iaVz*]. Then 

J d^xe,,ktriU-^d,UU-^d,UU-^dkU) (66) 

In the absence of the coupling to the gauge field, the Hopf 
term changes the spin and statistics of the skyrmion. 
Note that the gauge field couples to the skyrmion cur- 
rent density of the h vector. Therefore at long scales the 
skyrmion density is pinned at zero. Thus the resulting 
model is the three dimensional 0(3) vector model where 
skyrmion configurations have been suppressed. This 
leads to enlargement of the degrees of freedom to a four 
component field and the symmetry to 0(4). To see this 
most simply, rewrite the 0(3) sigma model above in the 
CP^ representation. We have 

^0 = j d'x-^\{d,-iAi)z\^ (67) 

Sj = / (fx-^CijkaidjAk (68) 
J 2.7: 

We may now integrate over the gauge field a. This is con- 
veniently done by choosing a gauge where Ai is transverse 
(ie V.A^ 0). We get 

S = S, + Sh (69) 

= f d^x-\ (a, - zA,) + (70) 

Thus the CP^ gauge field Ai has been rendered massive 
due to the coupling to the fluctuating gauge field a. At 
long distances and low energies we may thus drop the 
gauge field altogether. (Strictly speaking we should just 
integrate out the gauge field to generate a term that is 
quartic in the z fields and involves two derivatives - this 
term is expected to be irrelevant at the critical point of 
the resulting theory. A similar result was established in 
Ref.^^ in the absence of the topological term). The action 
Sz then becomes 

Sz= I (f'x-\diz\'^ (71) 
J 9 

which has 0(4) symmetry as can be made manifest by 
rewriting in terms of the real and imaginary parts of the 
two components of z. What happens of the Hopf term? 
Clearly Sh is invariant both under independent left and 
right multiplications by constant SU{2) matrices. Thus 
it too is 0(4) invariant. A simple calculation shows that 
the Hopf invariant H is exactly equal to Q where Q is the 
topological index characterizing configurations of an 0(4) 
unit vector field in three dimensions. Thus the action 
reduces to the 0(4) model at 9 = n. 



Note that in the 0(4) symmetry broken phase there 
are three gapless linear dispersing modes. In the gauge 
theory description this corresponds to broken chiral sym- 
metry with < n >^ 0. Here again there are three gapless 
linear dispersing modes - two are spin waves in n while 
the third is simply the gapless photon. 

A somewhat similar possible duality between massless 
N = 2 QED3 and the usual critical 0(4) model was 
conjectured recently in Ref. using very different argu- 
ments. Our derivation shows that such a duality indeed 
exists but necessarily includes the topological term. Does 
the topological term make any difference to the proper- 
ties of the model? We turn to this question in the fol- 
lowing subsection. 



A. Isotropic 0(4) model with a 9 term 

What may we say about the properties of the isotropic 
(i.e 0(4) symmetric) model at = tt from the analysis 
in this paper? First we expect that a phase with bro- 
ken 0(4) symmetry is stable for weak coupling. The 9 
term presumably plays very little role at low energies 
in this phase as it only affects topological configurations 
that cost large energy. At strong coupling paramagnetic 
phases where 0(4) symmetry is preserved are presum- 
ably possible. The arguments in previous sections show 
that a trivial featureless paramagnet cannot exist. To see 
this first assume that such a paramagnetic phase can in- 
deed exist. Then turn on a small anisotropy that breaks 
the 0(4) symmetry to 0(3) x Z2. Such an anisotropy 
will have no effect on the ground state of such a trivial 
paramagnet - but the arguments of Section^shows that 
such trivial paramagnetic ground states do not exist in 
the anisotropic model. Thus such states are forbidden 
for the isotropic model as well. What are the possibili- 
ties then? Gapped phases can exist if the ground state 
has topological order. Gapless paramagnetic states are 
also not precluded by these arguments. 

Similar and more interesting conclusions may also be 
drawn by considering weak 0(2) x 0(2) anisotropy (in- 
stead of 0(3) X Z2). First note that with this anisotropy 
the model has an extra Z2 symmetry that interchanges 
the two 0(2) symmetries. It is clear that this is 
exactly equivalent to the Motrunich-Vishwanath self- 
duality symmetrjii^ of the equivalent easy plane NCCP^ 
model. We can now discuss the phase diagram of the 
9 ~ TT 0(4) model with weak 0(2) x 0(2) anisotropy in 
terms of the known phases and phase transitions of mod- 
els with the same field content as the easy plane NCCP^ 
model. First consider the 0(4) ordered phase. With the 
0(2) X 0(2) anisotropy this becomes the location of a 
first order 'spin fiop' transition between phases that sep- 
arately break either of the two distinct 0(2) symmetries. 
What about phases that preserve 0(2) x 0(2) symmetry? 
Here there are two possibilities. First the self-dual second 
order deconfined critical line of the easy plane NCCP^ ap- 
pears as a paramagnetic phase of the sigma model. Thus 
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this is a concrete example of the possibility of a gapless 
strong coupling paramagnetic phase at 6* = tt in the 0(4) 
model (albeit with some 0{2) x 0(2) anisotropy). It is 
concievable that the self-dual critical manifold of the easy 
plane NCCP^ admits a special multicrtitical point with 
higher 0(4) symmetry. The existence of such a fixed 
point would imply the possibility of a gapless paramag- 
netic phase at 9 — TT with full 0(4) symmetry. The other 
possibility that preserves full 0(2) x 0(2) symmetry is a 
gapped Z2 topologically ordered paramagnet which core- 
sponds to the Z2 spin liquid allowed in the lattice easy 
plane spin- 1/2 antiferromagnet. This phase should per- 
sist with full 0(4) symmetry as well. 

In view of the above it seems clear that the transitions 
out of the ordered phases with 0(2) x 0(2) anisotropy 
would be in a different universality class from those in a 
model without the 6 term. With full 0(4) symmetry it 
thus seems rather likely that at = tt the transition is in 
a different universaity class from that at 6 = 0. 



V. DISCUSSION AND PROSPECTS 

In this paper we have provided a number of examples 
illustrating how the physics of two dimensional spin-1/2 
quantum magnets can in principle be described in 'su- 
perspin' sigma models that describe the slow competing 
orders. As such these are as close that one might get 
to a Landau-Ginzburg description of these competing or- 
ders. An important feature of these sigma models is the 
presence of topological terms which reflect the underly- 
ing quantum nature of the spins. From the conventional 
point of view it is these topological terms that com- 
plicate a direct application of naive Landau-Ginzburg- 
Wilson(LGW) thinking to quantum magnetism. Thus 
they are at the root of the failure of the LGW paradigm 
in describing various quantum phases and phase transi- 
tions. 

The sigma model formulations provide a potential al- 
ternate to the gauge theoretic descriptions that have been 
used thus far to describe the non-LGW physics. However 
as things stand we know even less about how to handle 
the effects of topological terms than we know about gauge 
theories. So it is at present not clear how useful the sigma 
model will be. 

Finally our results raise the possibility of such a sigma 
model description for stable algebraic spin liquids in two 
dimensions (of which the dRVB state popular in cuprate 
physics may be a possible example). These also have 
slow power law corelations in a number of physical ob- 
servables. The results on iV = 2 massless QED3 pro- 
vide some positive hints. Stronger evidence is the sigma 
model description of the deconfined Neel-VBS critical 
point. These may be thought of as a special kind of 
algebraic spin liquid that has one relevant perturbation. 
In the gauge theory description they have gapless bosonic 
matter fields coupled to a non-compact U{1) gauge field. 
As these theories seem to have sigma model descriptions 



perhaps their fermionic cousins do as well. Perhaps such 
descriptions might even be useful! 



Acknowledgements 

We thank M. Freedman, M. Hermele, O. Motrunich, 
and Diptiman Sen for useful discussions. This work was 
initiated at the Aspen Center for Physics Summer Prog- 
arm on "Gauge Theories in Condensed Matter Physics" . 
TS also acknowledges funding from the NEC Corpora- 
tion, the Alfred P. Sloan Foundation, and an award from 
the The Research Corporation. MPAF was supported by 
the National Science Foundation through grant DMR- 
0210790. 



APPENDIX A: DUALITY OF THE 9 = tt 
0{2) X 0(2) MODEL TO NCCP^ 

Proceeding as usual we solve the current conservation 
condition Eqn. 1531 bv writing 

J2 - V X A (Al) 

with A an integer living on the links < ij >. We have 
for the action 

S = Stl+SA + SaA+Sa (A2) 

Sa = uJ2i^ X (A3) 
p 

SaA = Y.a2.yxA (A4) 
IJ 

with Sti and Sa as before. This must be supplemented 
with the mutual statistics condition 

n^„(-i)''-=i (A5) 

p 

To handle ths we write 

A,j ^ 2A^ + (A6) 
with s = 0, 1 and A' an integer. Then we have 

Y[a,,{-ir^^l (A7) 
p 

This is solved by 

a^ji-iy^^ = a^aj (A8) 

with ai = ±1. The integer condition on A' can be im- 
plemented (softly) with a term 

- tccos {2ttA[j) = -tccos {7t{A.,j - s,j)) (A9) 
= —tcCTijaiajCos (irAij) (AlO) 
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We now separate out the longitudinal part of A: 

A = I^ + -V9c (AH) 

TT 

Collecting all the pieces of the action together and inte- 
grating out the gauge field 02 we see that is massive. 
We will therefore simply drop it from now on. We are 
then left with 

S = Sn+Stc + Sai (A12) 
Stc = ~tc (TjjCOS (V^c) (A13) 

Sal = i^^(vxai)' (A14) 



In writing the term Stc we have absorbed the ai into the 
^gj. The sum over the Z2 gauge field may now be 
straightforwardly performed. It generates a number of 
terms of which the most important have the structure 



-K [cos{'V6+ - ai) + cos(W_ + ai) ) (A15) 



where 6± = 6c ± (f>i- Together with the term Sai this is 
precisely the action for the easy plane NCCP^ theory. 
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